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New results on atmospheric neutrinos from Soudan 2 * 
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for the Soudan-2 Collaboration 

(Argonne National Laboratory, University of Minnesota, Tufts University, Western Washington, USA; 
Oxford University, Rutherford Appleton Laboratory, UK) 

Neutrino interactions recorded in a 5.1 fiducial kiloton-year exposure of the Soudan-2 iron tracking calorimeter 
are analyzed for effects of neutrino oscillations. Using contained single track and single shower events, we update 
our measurement of the atmospheric f M /i/ e ratio-of-ratios and find R = 0.68 ± 0.11 ± 0.06. Assuming this 
anomalously low R- value is the result of flavor disappearance via z/ M to v T oscillation, we select samples 
of charged current events which offer good resolution, event-by-event, for L/E v reconstruction. Oscillation- 
weighted Monte Carlo events are fitted to these data events using a \ 2 function summed over bins of log(L/E v ). 
The region allowed in the (sin 2 26, Am 2 ) plane at 90% CL is obtained using the Feldman-Cousins procedure: 
0.46 < sin 2 26» < 1.0 and 2.2 x 10" 4 < Am 2 < 2.2 x 10" 2 eV 2 . A small but relatively energetic sample of partially 
contained f M events has also been isolated. Their distribution in \og(L/E v i s ) relative to null oscillation Monte 
Carlo is compatible with z/ M to v r oscillation scenarios within the parameters region allowed by our contained 
events. 



1. DETECTOR; DATA EXPOSURE 

The Soudan-2 experiment is currently taking 
data using its fine-grained iron tracking calor- 
imeter of total mass 963 tons. This detector im- 
ages non-relativistic as well as relativistic charged 
particles produced in atmospheric neutrino reac- 
tions, ft is operating underground at a depth of 
2100 meters-water-equivalent on level 27 of the 
Soudan Mine State Park in northern Minnesota 
(northwest of Sudbury). The calorimeter's mod- 
ular design enabled data-taking to commence in 
April 1989 when the detector was one quarter of 
its full size; assembly of the detector was com- 
pleted during 1993. Data-taking has continued 
with 85% live time, even though dynamite blast- 
ing has been underway nearby for the MINOS 
cavern excavation since Summer 1999. The data 
exposure as of this Conference is 5.40 fiducial 
kiloton- years (kty). Results presented here are 
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based upon a 5.1 kty exposure. 

The tracking calorimeter operates as a slow- 
drift (0.6 cm//is) time projection chamber. Its 
tracking elements are meter-long plastic drift 
tubes which are placed into the corrugations of 
steel sheets. The sheets are stacked to form a 
tracking lattice of honeycomb geometry. A stack 
is packaged as a calorimeter module and the de- 
tector is assembled building-block fashion using 
these modules [jj]. The calorimeter is surrounded 
on all sides by a cavern-liner active shield array 
of two or three layers of proportional tubes . 

Topologies for contained events in Soudan 2 
include single track and single shower events 
(mostly za, and v e quasi-elastic reactions) and 
multiprong events. Flavor-tagging proceeds 
straightforwardly: An event having a leading, 
non-scattering track with ionization dE/dx com- 
patible with muon mass is a candidate charged 
current (CC) event of flavor; an event hav- 
ing a prompt, relatively energetic shower prong 
is a candidate v e CC event. Recoil protons 
of momenta greater than approx. 350 MeV/c 
are imaged by the calorimeter, allowing a much- 
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improved measurement of the incident neutrino 
direction, especially for sub-GeV quasi-elastic re- 
actions. 

2. ATMOSPHERIC v FLAVOR RATIO 

We measure the atmospheric neutrino v^/v e 
flavor ratio-of-ratios R using single track and 
single shower events which are fully contained 
within the calorimeter (all hits more than 20 cm 
from the nearest surface). These samples con- 
tain mostly quasi-elastic neutrino reactions, but 
include a background of photon and neutron reac- 
tions originating in cosmic ray muon interactions 
in the surrounding cavern rock. The latter "rock 
events" are mostly tagged by coincident hits in 
the active shield, however some are unaccompa- 
nied by shield hits and constitute a background. 
The amount of zero-shield-hit rock background 
in a neutrino event sample is estimated by fitting 
event vertex-depth distributions to a combination 
of tagged-rock and v Monte Carlo distributions. 
As expected, the fits show the background to be 
mostly confined to outer regions of the calorime- 
ter. Details of our analysis procedures for quasi- 
elastic events can be found in Refs. ||. 

The track and shower event samples for our 5.1 
kty exposure are summarized in Table |l|. Our full 
detector Monte Carlo simulation of atmospheric 
neutrino interactions is based on the 1996 Bartol 
flux for the Soudan site 8 . 



Table 1 

Soudan-2 track and shower event samples from 
the 5.1 fiducial kiloton-year exposure. 



Tracks 


Showers 


Data, raw 133 


193 


Monte Carlo events 1097 


1017 


(norm, to 5.1 kty) 193.1 


179.0 


Data, bkgrd subtr. 105.1±12.7 


142.3±13.9 



After correction for cosmic ray muon induced 
background, the number of single track events ob- 
served in data is less than the number of single 
shower events, whereas the null oscillation Monte 



Carlo predicts the relative rates to be other-way- 
round. Consequently the flavor ratio-of-ratios ob- 
tained is less than 1.0 and is anomalous: 

R = 0.68 ± O.ll(stat) ± Q.06(sys). 

This value is equal to the R value obtained last 
summer using 4.6 kty exposure 

3. SAMPLE FOR L/E MEASUREMENT 

The phenomenology for to v r oscillations is 
quite specific; neutrinos of muon flavor can meta- 
morphose and thereby "disappear" according to 
the equation 



\26) 



1.27 Aro 2 [eV 2 ] ■ L[km] 



E v [GeV] 



(1) 



Consequently it is optimal to analyze for neutrino 
oscillations using the variable L/E v . With the 
Soudan-2 calorimeter, measurement of event en- 
ergy for charged current reactions is straightfor- 
ward; we do this with resolution AE/E which is 
20% for i/„ CC's and 23% for v e CC's. To de- 
termine the neutrino path length L however, the 
zenith angle 6 Z of the incident neutrino must be 
reconstructed with accuracy. The path length can 
be calculated from the zenith angle according to 

L(9 Z ) = 

y/(R-d) 2 cos 2 6 z + {d + h) {2R-d + h) 
-(R-d)cos9 z (2) 

where R is the Earth's radius, d is the depth of the 
detector, and h is the mean neutrino production 
height. The latter is a function of v flavor, v 
energy, and 9 Z || . 

We select from our data an event sample suited 
to this measurement. We use a quasi-elastic track 
or shower event provided that the recoil proton 
is measured and that Pi ep t exceeds 150 MeV/c; 
otherwise, if the recoil nucleon is not visible, we 
require the single lepton to have E V i S great than 
600 MeV. We also select multiprong events, pro- 
vided they are energetic (E V i S greater than 700 
MeV) and have vector sum of P V i S exceeding 450 
MeV/c (to ensure clear directionality). Addition- 
ally, the final state lepton momenta are required 
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to exceed 250 MeV/c. For the selected sample, 
flavor tagging is estimated to be correct for more 
than 92% of events. The resolution for recovering 
the incident neutrino direction is evaluated us- 
ing the mean angular separation between "true" 
versus reconstructed neutrino direction in Monte 
Carlo events. The mean separations are 33.2° for 
2/ M CC's and 21.3° for v e CC's . The resolution 
in log LjE v {L in kilometers, E v in GeV) is bet- 
ter than 0.5 for the selected sample. Hereafter we 
refer to these events as "HiRes events" . 

The zero-shield-hit rock background, as esti- 
mated by the fits to event vertex depth distribu- 
tions, comprises 6.8% (5.1%) of the (v e ) flavor 
sample of HiRes events. 

Table |2| shows the HiRes event populations. Af- 
ter background subtraction there are 106.3 data 
events of flavor and 132.8 events of ^ e -flavor. 
Using these events, whose mean energy is higher 
than that of our track and shower events, the 
ratio-of-ratios is R = 0.67 ± 0.12, which is also 
significantly less than 1.0. 



Table 2 

Event samples selected for good LjE v resolu- 
tion, including atmospheric v data (without, with 
background subtraction) and v Monte Carlo sam- 
ples. The MC rates are shown normalized to the 
v e CC data. 









Data, raw 


114.0±10.7 


140.0±11.8 


Data, subt. 


106.3±14.7 


132.8±13.4 


Monte Carlo 


158. 5± 4.8 


132.8±4.4 



The atmospheric Monte Carlo (MC) sample 
represents 28.2 kiloton years of exposure. The 
MC event rates displayed in Table [| have been 
normalized to the v e data sample. This normal- 
ization is equivalent to a reduction of the Bar- 
tol neutrino fluxes by 21%. The assumption im- 
plicit with this adjustment is that the v e sample 
is devoid of oscillation effects. Figs. |], || and || 
show HiRes distributions with this normalization 
in place. 

Fig. [I] shows the distributions of these sam- 



ples in cosine of the zenith angle. For v e events, 
the shape of the distribution for data (Fig. [l]a,, 
crosses) coincides with that predicted by the 
Monte Carlo (dashed histogram) for null oscilla- 
tions. The distribution of data however, falls 
below the MC prediction in all bins (Fig. |l|b) 
with the relative dearth being more pronounced 
for v^s incident from below horizon. Distribu- 
tions in log (L/E„) for HiRes events are shown 
in Fig. ^. For null oscillations this variable dis- 
tributes according to a 'phase space' which re- 
flects the neutrino points-of-origin throughout the 
spherical shell volume of the Earth's atmosphere. 
That is, down-going j/'s populate the peak at 
lower \og(L / E„) from 0.0 to 2.0. Neutrinos inci- 
dent from/near the horizon occur within the dip 
region extending from 2.0 to 2.6, while upward- 
going neutrinos populate the peak at higher val- 
ues. Fig. ||a shows that, allowing for statistical 
fluctuations, the v e data follows the shape of the 
null oscillation MC distribution. In contrast, the 
data (Fig. ||b) falls below the null ocillation 
MC for all but the most vertically down-going 
flux. 

4. (sin 2 26>, Am 2 ) ALLOWED REGION 

To convert results of our atmospheric neutrino 
simulation generated under the no-oscillation hy- 
pothesis into simulated neutrino oscillation data, 
we apply to every MC event an L/E'y-dependent 
weight representing the probability of flavor 
survival for a given Am 2 and sin 2 29. 

An exploratory matchup of data with tri- 
al oscillation scenarios is shown in Fig. |3| For 
the mixing angle sin 2 29 fixed at unity, we plot 
the MC distribution weighted for to v T os- 
cillations with differing Am 2 settings. With 
Am 2 = 10~ 4 eV 2 (Fig. |a), the oscillation pre- 
diction lies above the data in most L/E v bins. 
With Am 2 = 10~ 3 eV 2 (Fig. |b), the upgoing 
v flux is now better described by the MC, al- 
though the expectation for horizontal and down- 
going neutrinos remains somewhat high. With 
Am 2 = 7 x 10~ 3 eV 2 , a rough agreement over- 
all is achieved (Fig. ||c). Going to higher Am 2 , 
we find that for Am 2 = 10 _1 eV 2 the oscillation- 
weighted MC falls below the data in most bins 
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Figure 1 . Distributions of cos 9 Z for v e and fla- 
vor HiRes samples. Data (crosses) are compared 
to the null oscillation Monte Carlo (dashed his- 
tograms) where the MC has been rate-normalized 
to the v e data. 



Figure 2. Distributions of log (L/E v ) for v e and 
charged current events compared to the atmo- 
spheric neutrino MC with no oscillations. The 
MC is shown rate-normalized to the v e data. 



(Fig. |3|d). These hints concerning the parameters 
regime preferred by the data are borne out by a 
more considered analysis, as we now show. 

To determine the neutrino oscillation param- 
eters Am 2 and sin 2 29 from our data, we con- 
struct a x 2 function over the plane-of-parameters. 
For points in the physical region of the 

(sin 2 29i, log Am 2 ) plane, we fit the MC expec- 
tion to our data at each point. The MC flux nor- 
malization, f v as well as sin 2 29 i and log Am 2 , is 
a free parameter: 

(X 2 data h= X 2 (sin 2 2^, Am 2 , f v ) 

A (N k (data - bkgd) - /„ • N k (MC)f 



We assume that the oscillation affecting our data 
is purely into v T and that the v e data is unaf- 
fected. 

The x 2 is summed over data bins containing 
our selected (HiRes) and v e samples, where 
k = 1 — 7 are log (L/E v ) bins, with k = 8 
containing all the v e events. The denominator cr 2 
accounts for finite statistics in the neutrino Monte 
Carlo and for uncertainty in the rock background 
in the v data. Not yet included are error terms 
which address systematic errors in the analysis, 
however preliminary examination shows statisti- 
cal errors to be the dominant error source in the 
analysis. The MC counts N k (MC) for the k th bin 
are constructed using oscillation weight factors. 

We find the location of minimum x 2 datal and 
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Figure 3. Comparison of log (L/E v ) distribu- 
tion for data (crosses) with expectations for 
Vp** v T oscillation with sin 2 29 = 1 (dashed his- 
tograms), for four different Am 2 values. 



P l0t ( A X 2 da tah Which is (XdatJij ~ (xLJmin- 

The A% 2 surface thereby obtained is shown in 
Fig. ^. A crater region of low \ 2 values is clearly 
discerned, at the bottom of which is a relatively 
fiat basin. The lowest point xtr-i* occurs at values 



sin 2 29 = 0.90, Am 2 = 7.9 x 10" 3 eV 2 , with flux 
normalization f v = 0.78. 

An additional structure is the A% 2 ridge which 
occurs at large mixing angle and for Am 2 above 
10 -2 eV 2 . For oscillation solutions in this regime, 
depletion in the downward-going v u neutrino flux 
with sub-GeV energies is predicted for to v T 
oscillations by equation (1) arising from the first 
oscillation minimum. Our HiRes events have suf- 
ficient resolution to show such an effect if it would 
be present. However, no pronounced depletion is 
observed, and so the \ 2 has a high value there. 

To find the region allowed for the oscilla- 
tion parameters by our data at 90% confidence 
level (CL), we use the method of Feldman and 




1 v&W^ 



Figure 4. The surface of A% 2 over the Am 2 , 
sin 2 29 plane; the MC normalization is allowed to 
adjust at each point. The oval at the bottom of 
the basin at large mixing angle denotes the Xmin 
location. 



Cousins 



At each of 2500 points 



(sin 29i, Amj) on a grid spanning the physical 
region of the plane parameters, we run 1000 simu- 
lated experiments. For each of the simulated sets, 
we find (Axgo)ij such that (Ax 2 im )ij is less than 
(Ax|o)ij f° r 90% of the simulated experiments at 
The surface defined by local Axg over the 
oscillation parameters plane is shown in Fig. |^. 
Note that the surface is not a plane at A^go = 
4.61, but rather has a concave shape. The central 
shaded portion is approximately A% 2 = 4.6, how- 
ever the outlying regions have A% 2 values which 
are lower. At each point over the physical region, 

if ( A X 2 datah is less than ( A xlo)i 3 , then be- 
longs to the allowed region of the 90% CL con- 
tour. 

The region allowed by our data at 90% CL 
is shown by the shaded area in Fig. |[ Al- 
though Xmin occurs at the location depicted by 
the solid circle, the relatively flat basin of our 
Ax 2 surface extends to lower Am 2 values. Su- 
perK has reported their best fit Am 2 value to be 
3.2 x f0~ 3 eV 2 ||; our data is compatible with 
that as well as with somewhat higher Am 2 val- 
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Figure 5. The surface of Ax| over the parame- 
ters plane described via the Feldman-Cousins pro- 
cedure. The intersection of this surface with the 
A\ 2 surface of Fig. [| defines the parameters re- 
gion allowed by Soudan-2 data at 90% CL. 



90% CL 

UZZj 68% CL 

(Feldman-Cousins) 



al 



sin 2 26 = 0.90 
Am 2 = 7.9x10 3 eV 2 
" (f v = 0.78) 

j i i I i i i I i i i L 




0.2 



0.4 



0.6 



_L 



0.8 1 

sin 2 (29) 



Figure 6. Soudan-2 allowed regions at 68% and 
90% CL with the best fit point for — > v T oscil- 
lations. 



ues. 

5. PARTIALLY CONTAINED EVENTS 

We plan to include more data in the above 
analysis. An additional data sample is comprised 
of flavor events which are partially contained. 
For each event of this category, the primary ver- 
tex is required to be > 80 cm (one hadronic in- 
teraction length) from exterior surfaces of the 
calorimeter, and the final state must contain a 
non-scattering, exiting track with ionization com- 
patible with a muon mass. This sample is useful 
because the assignment of flavor is reliable to 
better than 98%, and because the events are rel- 
atively energetic and consequently "point well" 
to the incident v direction. The mean energy for 
neutrinos which initiate PCEs is estimated to be 
4.7 GeV, to be compared to a mean energy of 
1.3 GeV for HiRes events. The mean angular 
deviation of the reconstructed direction versus 
the true direction (in Monte Carlo) is 14°. Unfor- 
tunately the number of PCEs is low, less than 
one-third the population of our HiRes sample. 



In order to isolate PCE two-prong and mul- 
tiprong topologies, the data events (with MC 
events interspersed throughout) are processed 
through a software filter; those which pass are 
scanned. The filter is designed to eliminate 
downward-stopping muons which have endpoint 
decays. In Soudan-2, an electron from a muon de- 
cay near rangeout can give rise to a small shower 
of < 10 hits from the end of a muon track; the 
topology is roughly akin to that of a neutrino- 
induced two-prong. Consequently care must be 
taken to avoid remnant up-down asymmetry in 
PCEs introduced by the filter. The problem is 
avoided by requiring that there be > 20 hits from 
a PCE vertex which are additional to the muon 
track. 

To mitigate against cosmic-ray induced back- 
grounds we require that any hit in the active 
shield which is coincident with a PCE, must be 
clearly associated with the exiting muon track. 
Occasionally it happens that a charged pion 
ejected from the cavern rock is incident upon the 
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calorimeter. If the pion penetrates by more than 
an interaction length and then scatters inelasti- 
cally, it can mimic a v PCE topology. Back- 
ground events of this type are removed by requir- 
ing the net momentum of the hadronic system of 
a PCE to lie within the same hemisphere which 
contains the candidate muon track. 
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Figure 7. Distribution of log (L/E V i S ) for par- 



tially contained events, 
events weighted for — ► 
sin 2 29 = 1 and with four different Am 2 values. 



compared to MC 
V oscillation, with 



With these selections we obtain 31 events for 
which the cosmic ray induced background is less 
than one event. The PCE rate predicted by 
our Monte Carlo with no oscillations is 40 events. 
The distribution of PCE data in log (L/E vis ) can 
be compared, as done previously for contained 
HiRes events, to representative oscillation sce- 
narios having sin 29 = 1.0. In Fig. @ we observe 
that, compared to the data (crosses) the oscilla- 
tion prediction (dashed histogram) is relatively 
high for Am 2 = 10 -4 eV 2 (Fig. fjja). However 
the prediction for Am 2 = 10 _1 eV^ (Fig. 0d) is 



too low relative to the data. Figs. 0b, c suggest 
that the scenario preferred by the data lies in the 
regime between Am 2 = 1 to 7 x 10~ 3 eV 2 . 

6. PLANS 

In the near future, we will include the par- 
tially contained events into our \ 2 fit to L/E. 
Additionally, a sample of upward-stopping muon 
events initiated by neutrino reactions below the 
detector has been isolated and will be analyzed 
for oscillation effects. And of course we will con- 
tinue to accumulate and analyze new data. Our 
goal is to keep the Soudan-2 detector running and 
tuned for the change-of-beams incident at our un- 
derground site. The change will be from atmo- 
spheric v's to Fermilab v^s in Fall 2003, at which 
time the detector will serve the MINOS experi- 
ment [g. 
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